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The object of this thesis was to investigate the stability 
of a ship with automatic steering in a following sea using the 
nonlinear equations of motion and to compare these results with 
the results of linear approximations. 

The solution of the nonlinear equations was obtained by a 
technique of treating nonlinear terms as linear terms with non- 
constant coefficients. Results showed that when drift angles 
become 0.1 radians or greater the nonlinear solution differs 
substantially from the linear solution. But the nonlinear solu- 
tion has a greater range of stability in the wave profile than the 
linear solution. This leads to the conclusion than a system 
designed using linear theory provides sufficient stability for 
the ship in a following sea. 

The most important effect which this approach neglected 
was the coupling effect of the roll or heel motion on the yaw 
motion. Before a complete and accurate solution to the steering 
problem can be found the coupling effect must be investigated and 
some method of expressing it analytically must be produced. 
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I NT RO DUCT I ON 



Steering a ship in a following sea and the danger of the 
ship broaching in a following sea are problems which have been 
prominent for many years. Before the invention of the autopilot 
the only solution was a skilled helmsman and an experienced 
captain who avoided the situation whenever possible. 

Steering in a following sea is not difficult until the speed 
of the ship becomes almost that of the waves . Then the forces and 
moments created by the waves on the ship are able to act for rel- 
atively long periods. In todays fast ships it is quite possible to 
encounter this situation of low frequency of encounter of ship to 
wave crest. There are two main solutions to the problem. The most 
obvious and by far the simplest to obtain is to reduce the speed 
of the ship which increases the frequency of encounter and short- 
ens the period during which the forces act. But in todays world, 
time is money and to slow the ship in following seas is undesir- 
able. 

With an autopilot the steering problem can be treated ana- 
lytically because the motion of the rudder is a function of the 
yaw and yaw velocity of the ship. The first explanations of the 
steering stability were made by using the linear terms of the 
equations of motion. However, in ship motions such as turning 
maneuvers the linear equations did not give adequate analyti- 
cal results. Research has been done with the nonlineal: equations 
of motion and solutions have been obtained using important non- 
linear terms of the equations .These solutions were found to give 
adequate predictions of ship motion maneuvers. 

In a following sea with small frequency of encounter the 
motions may become large and the linear solution may be very 
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inaccurate. Until now no one has tried to use the nonlinear 
equations to predict the stability of a ship traveling in a 
following sea. 

The solution of the nonlinear equations of motion in the 
case of following seas is possible only after making several 
simplifying assumptions. The velocity of the ship is assumed to 
be the same as the wave thus the frequency of encounter is zero. 
The ship position in the wave profile is constant. The forces 
and moments created by the waves on the ship are determined by 
the Froude-Krylov Hypothesis. Cross coupling terms of roll and 
pitch motion are neglected in the equations. The surging motion 
of the ship is neglected. Having made these four assumptions it 
is now possible to solve the nonlinear equations of motion for a 
ship in a following sea. 
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PROCEDURE 



NOMENCLATURE 

The motions of the ship are based on axes through the center 
of gravity of the ship,Figure I, and the terms are defined below 
or in the reference material. 

et = X-' V Heading angle from direction of advance of the waves 
Angle of drift or side slip 

Coefficient of yaw feedback in rudder control system 
Wave length 

Yaw angle or heading error 

Coefficient of yaw velocity feedback in rudder control 
system 

Horizontal distance between wave crest and the center of 
gravity of the ship 
Block Coefficient 

Lateral force exerted by the waves on the ship 
Yawing moment exerted by waves on the ship 
Wave amplitude 
Rudder angle 

The derivatives and the nondimensionalizing factors for 
them are as defined in references (6) , (10) and (11). 



B = 
* = 
X = 

r = 

<r _ 
a. = 

C k = 

y 

n = 
h 

S = 
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Figure I 

Definitions of Position, Distances 
and Angles 




* 
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NONLINEAR EQUATIONS OF MOTION 



The motions of the ship may be represented by nonlinear as 
derived in reference (10) . These equations have been nondimension- 
alized in accordance with reference (11) . The coefficients of the 
terms in the equations were determined for the Mariner hull model 
in reference (6) . Using these equations and including only terms 
which reference (6) determined were important or of minor impor- 
tance, the Y and N equations are expressed as 

Y Equation 

(m - Yy) v + (mx G - Y<f)y - Y v v -1/6 Y vvv v 3 - (Y y - mu)*- 
-1/2 Yy VV YV 2 - Y -1/6 Y sfS * 3 -1/2 Nj vv *v 2 = 0 

N Equation 

(mx G - N • ) v + (I z - Nt)y - N v v -1/6 N vvv v3 -(Nr~ 

(2) 

mx G u) 'f' “1/2 N y vv yv 2 - -1/6 N ssf 6 3 -1/2 N^ vv &v 2 

= 0 

Substituting £=-v into equations (1) and (2) equivalent 

equations are obtained using drift angle in lieu of lateral 

/ 

velocity. Adding to equations (1) and (2) the force and moment 
excitation of the waves as developed in Appendix B 

Y Equation 

(m - Y^) B + (mx G - Y y ) T - Y v 6 -1/6 Y vvv # 3 - (Y y - mu) *f> 
-1/2 Y^ vy *G 2 - Y 6 S -1/6 Y ifi $3 -1/2 Y* vv S 02 = Fyy5 ( 3 ) 
cos 2 > sin(2 ffa/ji ) - 'riFyy/i# 5 cos 2 -y sin(2 s’a/j. ) 

N Equation 

(mx^ - Ny) O + (I z -N ) V” “ N v G -1/6 N vvv <? 3 -(Nv~ 
fu-Cqu)*' -1/2 N j, vv V^ 2 - N s 6 -1/6 N fS j $ 3 -1/2 N s vv * e2 

(4) 

= 1/2 F xx 5 cos 2 <^ cos(2 77 a/?, ) - V 1/2 0F XX /^^ I cos 2 ^ cos ( 

2 7T a/x ) 

1 = k/Cb m’4 -rr 2 hL/* 2 (5) 
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For a ship with an autopilot having yaw and yaw velocity 
feedback control the rudder equation is 

<J+ t i = yy+ r4- (6) 

LINEAR SOLUTION OF EQUATIONS 

The three equations in three unknowns are sufficient to give 
a solution. The difficulty arises in the fact that the equations 
are nonlinear differential equations and an analytical solution is 
not known. As a first approximation to the solution neglect all 
nonlinear terms and find the linear solution. This results in the 
following equations 

Y Equation 

(m - + (mx G - Yf)f - (Y Y - mu)f* - Y v e - Y* S + 

, - ( 7 ) 

JFyy/^ i cos'V sin ( 2 7T a/ A ) = F yy I cos'V sin(2fi"a/>> ) 

N Equation 

(mx G - Ny) e + (I 2 - N Y ) *> -(N^- mx G u) v -N v G - N ^ i + 

1/2 dF xx /i-sp* cos 2 -* cos ( 2 *r a/^ )y = l/ 2 F xx 5 cos 2 ^ (8) 

cos (2 rra/^. ) 

Rudder Equation 

s + t i = y *r + r* ± O) 

Using letter designations for the coefficients and the values 
for the coefficients as determined in reference (6) the equations 
become 

AO+Cy + Dtf + E Y + G $ + II Y = I 

JO + K v + L <3 + My + P j + R Y = T (10) 

<£ + t $ = 

where the values of letter coefficients are given in Table I.. 
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Table I 



Values of Coefficients for Mariner Hull 



A = 


(m - Y^) 


= .01546 


J = (mx G - N^) 


= .000227 


C = 


(mx G - Y r ) 


= -.000086 


K = (I 2 - N*>) 


= .00083 


D = 


“Y v 


= .011604 


L = -N v 


= .00264 


E = 


- (Y y - mu) 


= .00499 


M = - (N <i> - mx G u) 


= .00166 


G = 


-Y* 


= -.002779 


P = -N<f 


■ .00138 



H = cos 2 ^ sin(2/ra/^) R = *5 iF xx /j* i cos 2 -y cos ( 2«"a/>. ) 

I = Fyy 5 cos^sin (2*a/>. ) T = % F xx 5 cos^ cos (2 to/* ) 

The general solution of the simlutaneous equations 

(As + D) &c + (Cs 2 + Es + H) Vi + G S„ = 0 
(Js + L) c L + (Ks 2 + Ms + R) r . + P Si =. 0 (11) 

- ( c- s + * ) Vi +(ts + 1) Si =0 

where = c, e s '*- + e t e s £ + oe. + o f e s «b 

Vi = ^e st + 

£ = & e s £ + 



can be found by using determiants. The rrots of the solution are 
thus the stability roots of the system. The roots can be found by 
finding values of s for which the determinant is equal to zero. 



CONTROLS FIXED STABILITY IN WAVES 



For controls fixed stability there is no rudder term because 
the rudder deflection is zero. The controls fixed stability is 
determined from the following equation: 

(AK - JC) s 3 + (DK + AM - LC - JE) s 2 



( 12 ) 



+ (DM + AR - EL - JH) S + (DR - LH) = 0 



7 



Substituting real values of coefficients from Table I into 
equation (12) the stability of the system is determined by use of 
Routh's Criterion for various positions of the ship in the wave 
profile and for various wave heights and wave lengths. Results are 
shown in Figures II and III. 

STABILITY WITH AN AUTOPILOT IN WAVES 

For Stability with an autopilot in waves the solution 
determined from the following equation: 

(AK - JC)t s 4 + ((KD + AM - JE - LC)t + (AK - JC) ) s 3 
+ ((DM + RA - JH - LE) t + (KD + AM - JE - LC) + (AP - 
JG) <7- ) s 2 + ( (RD - LH)t + (DM + RA -JH - LE) + (PD - 
GL ) r + (AP - JG) ¥ ) s + (RD - LH + (PD - GL) i* ) =0 
Substituting real values of coefficients from Table I 
equation (13) and assigning values of X + *~ feedback coefficients 
the stability of the system can be found using Routh's Criterion. 
Results are shown in Figures IV, V and VI for the stability of 
the system in the wave profile for varying values of yaw feedback 
coefficient, with various wave heights and wave lengths. 

In computing the data for Figures IV, V and VI it became 
obvious that yaw velocity feedback had no effect on the range of 
stability in the wave profile. Some method of showing the effect 
of yaw velocity feedback on the magnitude of stability was need- 
ed. This was done by assuming negligible time lag in the system. 
This makes one root very large and negative thus effectively re- 
ducing the equation to a cubic for which the solution of the 
important stability roots was obtained. In solving the cubic only 
the value of the largest stability root is determined. Reducing 
equation (13) to a cubic and solving, the effect of yaw velocity 
feedback on stability is as shown in Figure VII. 



is 



. (13) 



into 
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NONLINEAR TERMS 



Thus far the affect of yaw feedback coefficient and yaw 
velocity feedback coefficient has been shown using the linear 
equations of motion. Now some method of analytically expressing 
nonlinear terms must be developed in order to show their effect 
on the stability solution of the problem. 

Going back to equations (3) and (4), the nonlinear terms of 
any importance are found to be the £ 2 , B 3 and <5 & 2 terms. 
Treating the terms as linear with nonconstant coefficients they 
can be expressed as follows: 

Y vvv c 3 becomes ( e 2 Y vvv ) b 
y £vv becomes { e 2 Y,j vv ) $ 

etc. 

All the nonlinear terms can be treated in the above form. 
Putting these terms into equation (13) after assigning letter 
coefficients as shown in Table II the following determinental 
equation is produced: 

(AK - JC)t s 4 + ( (KD + KD ' O 2 + AM + AQ & 2 - JE 
- JF G > 2 - LC - L'C C 2 )t + (AK - JC) ) S 3 + ((DM + D'Mfl 2 
+ DQ C 2 + D'Q /3 4 + RA - JH - LE - L'E B 2 - LF B 2 

-L'F G 4 ) t + (KD + KD' B 2 + AM + AQ f? 2 - JE - JFB 

-LC-L'C£ 2 ) + (AP + AP , <? 2 -JG-JG'<? 2 )e')s 2 (14) 

+ ((RD + RD' B 2 - LH - L'H B 2 ) t + (DM + D'M 3 2 + DQ C 2 

+ DQ £3 4 + RA - JH - LE-L'E © 2 - LF B 2 - L'F C 4 ) + (AP 

+ AP' 0 2 - JG - JG' G 2 ) * + (PD + P'D 3 2 + PD' G 2 + 

P ' D ' 8 4 - GL - GL ' 3 2 - G ' L <3 2 - G'L'G 4 )^) s +((RD 

+ RD' 3 2 - LH - L'Hfi 2 ) + (PD + P ' D G 2 + PD' 3 2 + 

P'D' G 4 - GL - GL % G> 2 - G'L8 2 - G ' L ' <3 4 ) ^ ) = 0 
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Table II 



Coefficients of Nonlinear Terms for 



Mariner Hull 



F — —1/2 Y yy 


= 


-.0765 


Q = -1/2 N ^ vv 


= .0274 


D’= -Y vvv 


= 


.0808 


L 1 = 


= -.0164 


G ’ = - Y <SW 


= 


.0119 


P *= - N 6W 


= .00489 


STABILITY WITH AUTOPILOT 


IN 


WAVES INCLUDING NONLINEAR TERMS 


Substituting values 


from Tables 


I and II into 


equation (14) 



and by choosing an instantaneous value of drift angle, Q, a 
solution for the stability can be determined for that instant of 
time. By choosing a value of drift angle which is felt to be the 
largest expected value, the solution for the roots is obtained. 
Plotting this solution along with the linear solution a range in 
which the true solution lies is found. These results are plotted 
in Figures VIII and IX. 

From all of the previous equations it is now possible to 
investigate the steering stability of any ship in waves and to 
determine what steering characteristics are necessary to have 
directional stability in a following sea at any given sea con- 
dition. 

STABILITY FOR MARINER HULL IN REAL SEA CONDITIONS 

In order to obtain information necessary to design a steer- 
ing system that will ensure directional stability at all times the 
effects of time lag, feedback coefficients and rudder area on the 
stability must be known. There may be limitations on the values of 
yaw feedback coefficient and yaw velocity feedback coefficient 
because operation in head seas and calm water must be part of the 
design. These operations may require some range of values of y and 
Time lag does not appear to be a critical parameter. Present 
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system time lags are in the area of 0.1 and this is sufficiently 
small for stability in following seas. 

The first consideration is that the rudder area must be suf- 
ficient to counteract the moment created by the waves. For the 
Mariner hull a check for sufficient rudder area is made by equa- 
ting the moment created by the waves to the moment which can be 
applied by the rudder. Figure X shows the angle needed to over- 
come moment of waves alone. The present Mariner rudder appears to 
be large enough to counteract wave moment. 

In all of the developments using the Mariner hull the assump- 
tion is made that the ship is somehow able to maintain the same 
speed as the waves. Thus the problem as developed represents a 
more dangerous condition than is actually present. 

The next step is to find the value of yaw angle feedback 
coefficient necessary to make the system stable in design sea 
conditions. Using equation (13) and the values of Tables I and II 
the feedback coefficient necessary for stability is found. Re- 
sults are shown in Figure XI. 

Another question that arises is how would changing the rud- 
der area change the stability feedback coefficients. This is done 
by using equation (13) and varying the rudder force and moment 
derivatives. Using the rudder coefficient changes of reference (9) 
as a basis a direct percentage change in Mariner coefficients is 
made in order to obtain an estimate of how Mariner coefficients 
would change with the same rudder area changes. The effect of 
rudder area changes on angle to overcome wave moment are shown in 
Figure XII. The effects of rudder area changes on yaw angle feed- 
back coefficient are shown in Figure XIII. 

From the results of these investigations in following seas 
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and from investigations in head seas and calm water maneuvers th 
best combination of rudder area, time lag, and steering system 
feedback coefficients can be chosen. 
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RESULTS 



The following figures show how the directional stability of 
a ship is affected by the presence of waves when running in a 
following sea. 

Figures II and III show how the range of instability in the 
wave profile for the ship with no steering control for varying 
wave height and wave length. 

Figures IV, V and VI show how the range of instability can be 
decreased and even eliminated by the addition of automatic steer- 
ing controls with ya w angle feedback control. 

Figure VII shows how the real and imaginary part of the 
largest stability root varies in the wave profile for a system 
which is stable all the time for two values of yaw velocity feed- 
back coefficient. 

Figures VIII and IX show how the nonlinear terms in the 
equations of motion affect the largest stability root and the 
range of stability in the wave profile. 

Figure X shows the rudder angle necessary to counteract the 
constant wave moment acting on the Mariner hull. 

Figure XI shows the minimum required yaw angle feedback 
coefficient ncessary for a stable system as a function of the 
sea condition. 

Figure XII and XIII show how the rudder angle reuired to 
counteract wave moment and the required yaw angle feedback coef- 
ficient are changed when rudder area is changed. 

Figure XIV shows how the minimum required yaw angle feedback 
coefficient varies as the size of the ship changes for operation 
in the same sea conditions. 
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Figure II 



Stability of Controls Fixed Ship in VJaves 




Figure III 

Stability of Controls Fixed Ship in Waves 
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/AW ANGLE FEEDBACK COEFF. <» WWANCLf FEEDBACK Coe/^F. ( V) 



Figure IV 



Stability of Autopiloted Ship with Yaw Angle 
Feedback Control 




Figure V 



Stability of Autopiloted Ship with Yaw Angle 
Feedback Control 
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Figure VI 



Stability of Autopiloted Ship v/ith Yaw Angle 
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P/JrtT RE/1L P/JBT 



Figure VII 

Real and Imaginary Parts of Largest Stability 

Root 

VL = 2 h/X = 1/40 
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YAW ANCLE FEEDBACK COtFF.O)) REAL PART 



Figure VIII 



Real Part of the Largest Stability Root 
<T= 1 1 >/L = 2 h/A = 1/40 & = .1 




Figure IX 

Effect of Nonlinear Terms on Range of Stability 

<r= o t = .1 >/L = 1 Va = 1/40 0 = 0.1 
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Figure X 



Rudder Angle Needed to Overcome Wave Moment 
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Figure XI 

Yaw Angle Feedback Coefficient Necessary for 
for Stability Using Linear Theory 
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Figure XII 

Rudder Angle Required to Counteract Wave Moment 
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Figure XIV 



Yaw Angle Feedback Coefficient Required for 
Stability for Various Ship Lengths Using 






Linear Theory 
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DISCUSSION 0? RESULTS 



NONLINEAR EQUATIONS 

The solution of the nonlinear equations presented an 
impossible task. To date the only solutions to the equations have 
been by numerical analysis which was designed to predict the 
motion of the ship during standard maneuvers in calm water. To 
obtain a numerical solution for steering in following seas was 
entirely beyond my comprehension. 

The only path left open was to try and find some method of 
obtaining an approximate analytical solution. Previous analytical 
solutions had been obtained using only the linear terms of the 
equations. Investigation of the important nonlinear terms resulted 
in the discovery that they were all a function of drift angle 
squared or drift angle cubed. By treating the drift angle squared 
as a non constant coefficient equations evolved as linear with 
non constant coefficients. These equations were solved by choos- 
ing a maximum expected drift angle and putting it into the equa- 
tions as a constant. This gave a stability solution for an 
instantaneous value of drift angle. The maximum expected drift 
angle solution is then the solution for largest ship motions. 

Comparing the linear solution and nonlinear solution of the 
largest stability root,Figure VIII, the nonlinear solution shows 
less stability than the linear solution. Figure IX shows that 
the ship has a greater range of stability in the wave profile 
than predicted by the linear solution. These results correlate 
with the nonlinear solutions of the standard maneuvers. In the 
standard maneuvers of reference {9) the linear solution predicts 
less stability and thus better turning characteristics than is 
actually present. The nonlinear terras represent forces and moments 
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which are always opposing the yawing motion. As a result the 
range of stability in the v/ave profile is actually greater than 
the linear prediction but the oscillating motion of the system 
does not decay as rapidly as linear theory predicts. 

The actual stability of the system lies somewhere in 
between the linear prediction and the nonlinear prediction 
because the drift angle is always changing as the ship moves 
through the waves. The graphs are plotted for a drift angle of 
six degrees. In the turning circle maneuver the maximum drift 
angle was nine degrees and for the zig zag maneuver the drift 
angle reached seven degrees. 

In designing a steering system the use of the linear theory 
will provide more than adequate control and so the design is 
sufficient for all conditions. 

DESIGN OF STEERING CONTROLS 

In the design of an automatic steering system the rudder 
size, steering equipment size and the feedback to the control 
circuits are the parameters which must be determined. It is very 
important to be able to see what effect these parameters have on 
the steering problem in a following sea. 

The rudder size and the steering gear size are initially 
determined by the requirements set for the ship to make the stan- 
dard maneuvers in calm water. 

With this rudder size Figures II and III show the stability 
region for the ship if it has no controls and the rudder remains 
fixed at zero deflection. The ship is unstable on the negative 
slope of the wave and stable on the positive slope. This stability 
region is unchanged by v/ave height but does shift with varying 
wave length because the v/ave slope changes. 
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Adding automatic controls to the ship with its present 
rudder - . Figures IV, V and VI shov; how the value of yaw angle 
feedback coefficient required for stability varies for different 
types of following seas. The value of yaw angle feedback coeffi- 
cient needed for stability is very dependent upon the size and 
length of the waves. As wave height increases the more yaw angle 
feedback is needed to make the system stable. The important wave 
lengths are in the tange 1.0 to 8.0 of ship length. Shorter than 
this and the wave has little affect on steering. Longer waves 
require the same or less feedback due to decreasing slope of the 
wave . 

Yaw angle feedback is the controling feedback for stability. 
The effect of adding ya w velocity feedback to the system is shown 
in Figure VII. At zero frequency of encounter with the waves it 
has an unstabilizing affect. But this is for our own hypothetical 
situation. When the frequency of encounter is not zero but has a 
value of 0.5 radians/sec. or larger then L.J. Rydill , reference (8) 
has shown that velocity feedback helps stability. So in an actual 
following sea situation a value of yaw angle velocity feedback is 
important for good ship control. 

Taking the Mariner hull in a following sea situation the 
effect of wave forces on stability can be shown. The effect of 
various sea states on a 500 foot Mariner is shown in Figures X 
and XI. Figure XI shov/s the minimum value of yaw angle feedback 
coefficient for operation in different sea states. Figures XII 
and XIII show how the minimum value of yaw angle feedback coef- 
ficient is changed by changing the rudder area. Figure XIV show 
how the required yaw angle feedback coefficient differs for a 
200 foot and a 100 foot Mariner as compared to the 500 foot ship. 
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Having obtained the data for required yaw angle feedback 
coefficient with rudder area variation, this must be added to data 
generated for steering in head seas and in standard maneuvers to 
give a true picture of the limitations of the steering parameters 
in the different modes of operation. With this data one should be 
able to design the best possible steering system. 

The results for the real ship case are highly dependent on 
available wave spectrum data. The data used here was for fully 
developed ocean waves as defined by the U.S. Hydrographic Office. 

In coastal waters it is very possible that the wave height to wave 
length ratio will be much larger for given wave lengths due to 
shallow water effects on the waves. This fact would require much 
larger values of yaw angle feedback coefficient to stabilize the 
ship. In designing a steering system it is most important that 
accurate data on operating sea conditions be used. 

LIMITATIONS OF DESIGN PROCEDURE 

There are several conditions and effects which may be of 
importance, that are not accounted for in this investigation. The 
most important is the cross coupling effect of the roll motion on 
the yawing motion. No one has yet been able to express analytical- 
ly the effect of rolling motions on the steering problem. It is 
easy to see that when the ship rolls or heels during yawing motion 
that it becomes unsymetrical about the reference centerline and 
side forces and yawing moments are created. The larger the roll 
or heel angle the greater the forces and moments. This effect is 
greatly dependent on the static stability of the transverse motion. 
Investigation by the Germans, reference (12) ,in capsizing model 
tests shows that as the static stability of the ship is reduced 
it becomes less stable in yawing motions until it reaches a point 
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where the model broaches and eventually capsizes. From their 
investigation the importance of the yaw- roll coupling is shown. 
Until someone develops a method to express this effect a complete 
solution to the steering problem is not really known. 

A condition which has not been investigated is the condition 
where the frequency of encounter with the waves is such that the 
wave moments are exerted when the ship is yawing away from the 
prescribed course such that the motion appears as a spring mass 
system being excited at half frequencies of its natural frequency. 
It seems possible that the ship yawing motion may somehow be 
expressed as a simple spring-mass-damper system with some natural 
frequency. Then the system could be checked at the half frequen- 
cies to see if the motion is stable or not. By changing the 
ammount of feedback the natural frequency would be changed and 
data could be obtained which may place further limitations on 
the values of the feedback coefficients. 

Another effect which has been neglected is that of the 
wind force and moment on the ship. There may, in a sufficiently 
strong wind, be forces large enough to affect the steering 
problem. 

The effect of the surging motion of the ship in a following 
sea is another condition which may have important adverse affects 
on the steering problem. For this investigation the surging 
motion was neglected in order to simplify the mathmatics of the 
problem. The surging may cause sufficient changes in forces and 
moments so as to cause instability. 
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CONCLUSIONS 



1. The nonlinear equations of motion will provide a more 
accurate solution when operating in a following sea because drift 
angles can be expected to be large. 

2. The nonlinear solutions show that when designing a steering 
system use of the linear analysis will provide more stability than 
is actually required for operation in a following sea. 

3. At very low frequencies of encounter the yaw angle feedback 
in the only control needed for stability. Yaw velocity feedback 
at low frequencies has an unstabilizing effect but is important 
for good stability when the frequency of encounter is above 

0.5 radians/ sec.. 

4. Values of feedback coefficients can be determined in order 
to make the system stable when operating in following seas. 

5. Values of feedback coefficient are highly dependent on the 
wave spectrum in which the ship is designed to operate. 

6. The roll motion influence on the yaw motion can be very 
large for ships with low tranverse stability. 

7. The point where wave excitation is at half frequency of yaw 
motion natural frequency may be as dangerous a position as the 
zero frequency of encounter point. 



27 



RECOMMEN DAT I ON S 



There are three areas of the problem which need further 
investigation. 

1. Explanation of the Roll-Yaw coupling effect which will 
produce the greatest changes in the yaw motion problem. 

2. Investigation of half frequency excitation which may also 
produce some startling results. 

3. Investigation of the surging motion effects which may cause 
slight changes in the problem. 

It may also be possible to somehow produce the ship motions 
using an analog computer. If this can be done then the nonlinear 
solutions can be produced exactly and the enswer to the half 
frequency problem could be solved at the same time. 
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APPENDIX B 



LATERAL FORCE AND MOMENT EXCITATION FROM WAVES 

The moments and forces acting on the ship due to the waves 
are based on the Froude-Krylov Hypothesis which states that: 

a. The presence of the ship does not disturb the pressure 
distribution in the waves. 

b. The pressure distribution on the ship is determined by 
the hydrostatic pressure due to the position relative 
to the wave crest with a correction in order to take 
into account the decreasing orbital diameters of the 
water particles with increasing water depth. 

Assuming a symetric ship with vertical frames at the water- 
line and following the development in reference (1) the lateral 
force, y, and the yawing moment, n, due to wave pressure are: 



y = Fyy • LBT • / 3 g 2 tj h/> sin(2rr a />. - w e ) 

Having assumed the frequency of encounter is zero. w e = 0 
y = Fyy « LBT * /=>g 2?? h/^ sin(2 7/a/\ ) 
n = F xx * %L * LBT * x>g 2 n h//> cos (2^ aA ) 

The coefficients F yy and F xx are:. 



F = 2/ LB sin 'X £f w (x) cos(2'‘x/a ) dx 



- v*. 



F xx = 4/L^B sin rye ^f w (x) sin(2^x/A ) dx 

Because is not the true angle between wave direction and 
the centerline plane of the ship the equations must be modified as 
follows : 



True angle = ^y~ V 

y = (Fyy + cFyy) LBT * p g 2 R/a sin(2^a/A ) 
n = (F xx + <sF n ) hL • LBT * eg 2 w h/A cos (2 tt a/> ) 
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* F yy = -v 

A ^xx = ” 't' ^F xx / J rf. 

y = (Fyy - v iFyy/}^ ) LBT • /»g 2 w hA sin(2n-a/x ) 

ft = (F xx -f^F x X /j-y ^ * LBT * 2 tt hA cos( 2 t a/A ) 

Non dimensionalizing the equations, 
y* = y/h u 2 S 
n’ = HA U 2 SL 

Assuming that the ship velocity in the direction of the 
waves is the same as the wave velocity 
wave speed c = VgA/2 v 
ship speed u = c/ cos -y 

u 2 = gA / 2 tt cos 2 ^ 

y' = (Fyy - JFyy/i^ ) I cos 2 /> sin(2 7T a/A ) ' 
ft' “ ( F xx " V iF xx A^ ) h I cos 2 /y cos (2 tt aA ) 

I = k/C^ m' 4 T7 2 hL/A 2 

k is the factor used to take into account the influence of the 
decreasing orbital velocities as shown in reference (1) . 

On the following page are the graphs showing the values of 
Fyy and F X x coefficients for the Mariner hull. 
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Figure XV 



Coefficient for Sway Exciting Force 




Figure XVI 

Coefficient for Yaw Exciting Moment 
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Stability of autopiloted ships in a foil 
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